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Abstract. Dynamical contraction of a slowly-rotating magnetized cloud 
is studied using the magnetohydro dynamical (MHD) simulations. In the 
isothermal stage (n <^ n\ ~ 10 10 cm -3 ), the cloud evolves similarly to 
that expected from Larson-Penston self-similar solution and experiences 
the run-away collapse. However, after the central density exceeds ~ n&, 
an accretion disk is formed around the adiabatic core. Just outside the 
core, outflow is ejected by the effect of magnetic torque (magnetocen- 
trifugal wind). Since ~ 10% of the mass is ejected with almost all the 
angular momentum, the specific angular momentum of the protostellar 
core reduces to that observed in main-sequence stars. 



1. Introduction 

Fragmentation of the clouds (cloud cores) is thought to be affected by the amount 
of angular momentum in it. On the angular momentum, there has been a prob- 
lem that the specific angular momentum of new born stars (e.g. T Tauri stars) 
~ 6 x lO 16 (i2 ! „/2i?0) 2 (i :, /lOday) _1 cm 2 s~ is much smaller than that of the 

parent molecular cloud cores j c \ ~ 5 x 10 21 (.R/0.1pc) 2 (f2/4kms _1 pc _1 )cm 2 s . 
Here we assumed the observed velocity gradient comes from the rotational mo- 
tion. Further, the specific orbital angular momentum of a binary whose separa- 
tion is about R ~ 100AU is equal to j orbit =s 4 x 10 19 (i?/100AU) 1 / 2 (M/lM o ) 1 / 2 

cm 2 s If the angular momentum is conserved in the star forming process, 
neither binary nor single stars are formed. 

Angular momentum is reduced by the effect of magnetic fields (B-fields). 
Considering an oblate spheroidal cloud is threaded by poloidal B-fields parallel 
to the symmetry axis, the rotation motion generates torsional Alfven waves in 
the inlercloud medium, which propagate away from the cloud. The torsional 
Alfven wave causes the intercloud medium to rotate. In this case, the angular 
momentum in the cloud decreases with a timescale as tn = — § — n a -. r = 

D Q 2p a Va 

2i *Bq (47rGp a )~ 1 / 2 , where fi, a, p a , Va, and Bq represent, respectively, the an- 
gular rotation speed, the column density of the cloud, the ambient density, the 
Alfven speed in the ambient medium, and the magnetic flux density outside the 
cloud (Ebert et al. 1960). Since the first factor 2irGo-/Bo is larger than unity 
for magnetically supercritical clouds, the timescale ts is comparable with the 
free-fall timescale in the low- density ambient medium and thus the angular mo- 
mentum transfer may seem inefficient as long as the cloud contracts dynamically. 

Molecular outflow begins at some epoch in the contraction stage of the 
cloud. The molecular outflow ejects matter and thus angular momentum. Con- 



1 



2 



Kohji Tomis&ka 



sider the situation that 90 % of the inflowing matter forms a star and the rest 
10 % is ejected by the molecular outflow. On the other hand, 0.1% of the an- 
gular momentum goes to the newly formed star and the rest is ejected by the 
outflow. If this occurs, the specific angular momentum of a newly formed star 
is reduced by a factor 10~ 3 from that of the parent cloud. To confirm whether 
this mechanism works in a process of the star formation, we have studied the 
evolution of a slowly rotating cloud and the formation of molecular outflows. 

2. Numerical Method & Model 

We assume the gas is well coupled with the B-fields. The basic equations to solve 
are the ideal MHD equations with cylindrical symmetry and the Poisson equation 
for the gravitational potential. Since the gravitational collapse proceeds heavily 
nonhomologously, we have to cover the wide dynamic ranges for both the density 
and the size. This was attained by using the nested grid method, in which a 
number of grids with different spatial resolution are prepared (here we used 24 
grids from the coarsest L0 to the finest L23) and calculated simultaneously. The 
finer grids are for small-scale structure near the core while the coarser ones cover 
the global structure. 

As the initial condition, we assume a cylindrical isothermal cloud in hydro- 
static balance. B-field strength is taken as proportional to p 1 ^ 2 and the cloud's 
rotation axis coincides with the axis of the cylinder. Density perturbations with 
A equal to that of the most unstable mode are added. 

3. Run-away Collapse Stage 

Isothermal clouds experience the run-away collapse, in which the central density 
increases very rapidly in a finite timescale. This is common for spherical clouds 
without rotation and B-fields (Larson 1969) and for a disk formed in a rotating 
isothermal cloud (Norman et al 1980). The gas falling down along the B-fields 
forms a contracting disk (pseudo-disk) even for rotating and magnetized clouds. 
Contraction continues if the equation of state remains isothermal. In this phase, 
there is no evidence for the outflow (Fig.la,c). 

4. Accretion Stage 

After the central density reaches pA ~ 10 10 cm~ 3 , the thermal radiation from 
dust grains is trapped in the central core and an adiabatic core is formed (Larson 
1969). Since the dynamical timescale of the core becomes much shorter than 
that of the outer isothermal envelope, the isothermal gas begins to accrete to 
the nearly hydrostatic adiabatic atomic core (first core) . To study the evolution 
of this accretion stage, we use a composite polytrope model as below pa the gas 
is assumed isothermal; for p > pa the gas becomes adiabatic and the polytropic 
index T is assumed as 5/3 (see Tomisaka 1998). 

Figure 1 illustrates the change from the isothermal run-away collapse (Figl.a, 
c) to the adiabatic accretion phase (Fig. lb, d). Figure la shows that B-fields 
(red lines) run almost perpendicularly to the disk (blue surfaces) in the run-away 
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collapse phase. Toroidal B-fields and rotational motion are not predominant 
(Fig.lc). On the other hand, in the accretion phase B-fields are squeezed by the 
effect of inflow motion in the disk and they form an hourglass-shape configura- 
tion. Further, toroidal B-fields are formed by the rotational motions which are 
amplified in the accretion phase (Fig. Id). In this phase, outflow occurs (Fig. lb, 
d) at the distance ~ 50 — 100 AU from the center. The time elapsed between 
these two is only about 1000 yr. 

The Lorentz force F = (c/47r)j x B contains the toroidal component = 
(c/47r)j P oi x Bp i when there exists the poloidal current (j po i) which makes the 
toroidal B-fields B^. Therefore, it is shown that the magnetic torque N = r F^, 
which changes the distribution of the angular momentum, works effectively only 
after the toroidal B-field is amplified by the rotation motion. The toroidal B- 
fields become predominant after the core formation epoch. Further, Blandford 
& Peyne (1982) showed that efficient angular momentum transfer from the disk 
to the (cold) wind material occurs only when the angle between the B-field line 
and the disk is smaller than 60° . This condition is satisfied only in the accretion 
phase (compare Fig. la and. lb) 



5. Angular Momentum Transfer 



To study the effect of angular momen- 
tum transfer by the magnetic torque, 
specific angular momentum j is plot- 
ted against the accumulated mass from 
the center M in Figure 2 (Tomisaka 
2000). The difference between the initial . 
stage (dashed line) and the core forma- 
tion epoch (dotted line) does not come 
from the angular momentum transfer 
but from the segregation of the mass 
with small j to form the high-density 
central part. True transfer occurs af- 
ter the core formation. The magnetic 
torque increases j of the matter on the 
disk surface and the matter with excess 
j is ejected as the outflow, while the Figure 2: Specific angular momentum j 
magnetic torque decreases j inside the is plotted against the accumulated mass 
disk and this leads to accelerate the ac- from the center M. Open circles, rilled 
cretion. squares, and triangles denote, respec- 

That is, B-fields changes the an- tively, the beginning and the end of the 
gular momentum distribution along one run-away collapse phase and the accretion 
field line. As a result, in ~7000 yr after phase. 

the core formation M ~ 0.2M Q is accumulated in the core and j of the core be- 
comes as small as ~ 10 16 cm 2 s 1 (solid line). Therefore, it is clearly shown that 
the angular momentum problem of newborn stars is solved taking into account 
the molecular outflow driven by the magnetic torque formed just outside the 
first core (R & 100AU). 
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Figure 1. The left panels illustrate the structure in the run-away 
phase, while the right ones are for the accretion phase. The upper 
panels represent nearly edge-on views and the lower ones are pole- 
on views. These represent the structure captured by L10 grid and 
the diameter of the numerical box (white circle) is equal to ~ 300 
AU. Isodensity surfaces and B-field lines are illustrated as well as the 
velocity vectors and density countours. 



